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Tidal dissipation in rotating fluid bodies: a simplified 
model 
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ABSTRACT 

We study the tidal forcing, propagation and dissipation of linear inertial waves in a 
rotating fluid body. The intentionally simplified model involves a perfectly rigid core 
surrounded by a deep ocean consisting of a homogeneous incompressible fluid. Cen- 
trifugal effects are neglected, but the Coriolis force is considered in full, and dissipation 
occurs through viscous or frictional forces. The dissipation rate exhibits a complicated 
dependence on the tidal frequency and generally increases with the size of the core. In 
certain intervals of frequency, efficient dissipation is found to occur even for very small 
values of the coefficient of viscosity or friction. We discuss the results with reference 
to wave attractors, critical latitudes and other features of the propagation of inertial 
waves within the fluid, and comment on their relevance for tidal dissipation in planets 
and stars. 
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Tidal interactions determine the fate of short-period extra- 
solar planets. They affect their spin and orbital parameters, 
and in extreme cases may lead to the destruction of planets 
as a result of orbital decay or intense tidal heating. Tidal 
evolution is also important in close binary stars and in the 
satellite systems of the planets of the solar system. 

The efficiency of these processes is often parametrized 
by the tidal quality factor Q, which is a dimensionless 
inverse measure of the dissipative prop erties of the body 
consi dered as a forced oscillator (e.g. iGoldreich fc Soterl 
1966). The tidal Q of fluid bodies such as stars and 
giant planets is difficult to calculate from first principles, 
even in linear theory. The tidal disturbance generally 
consists of two parts. One (the 'equilibrium tide') is a 
quasi-hydrostatic bulge that is carried around the body 
by a smooth velocity field, while the other (the 'dynam- 
ical tide') consists of internal waves that are excited by 
the low-frequency tidal forcing and may have a short 
wavelength. Dissipation of the equilibrium tide can occur 
through its interaction with turbulent convection, but the 
damping rate is uncertain, especially when th e tidal period 
is sho r t compared to the co n vectiv e timescale (IZahnl 1966bl. 



tides one should consider the excitation, propagation and 
dissipation of low-frequency internal waves in rotating, 
stratified fluids. Calculations have been made for a variety 
of object s, using different s i mplifi c ations and approxima- 
tions b y IZahnl iil970l, Il975l. Il977l). ISavoniie fc Papaloizoul 
ill983l. 1 19841). ISavoniie. Papaloizou fc Alberts! (1 19951). 
ISavoniie fe Papaloizoul (1993), iPapaloizou fc bavonuc 
' 19971) . iLubow. Tout fc Livid (Il997l ). iTerquem et afl 
1998 ! ) . IGoodman fc Dickson! ll 19981) . IWitte fc SayonTiel 
1999 ) ISavoniie fc ; Wittd (120021), lOgilvie fc Lin! (|2004 
20071 ). IWul (l2005allb h IPapaloizou fc Ivanovl (|2005l ) and 



Ivanov fc Papaloizoul l|2007l ). This approach has perhaps 



19771 ; IGoldreich fc Keelevl 
19771 ; IGoodman fc Oh! 1 199 



19771; IGoldreich fc Nicholson! 
Penev et all 120071 ). This 



been most successful for early-type stars, in which gravity 
(or inertia-gravity) waves are excited near the base of the 
radiative envelope and propagate outwards until they are 
dissipated by radiative damping; in this case estimates can 
be made rather simply for the dissipation rate (or tidal 
torque) that are ind ependent of the details of the wave 
damping mechanism (|Goldreich fc Nicholsoniri989l ). 

In recent wor k on dynamical tid e s in rotating giant 
planets and stars (|Ogilvie fc Linll2004 I2007D we have em- 
phasized the role of inertial waves in convective regions, 
which are nearly adiabaticall y strat ifie d and do not support 
gravi ty waves (see also IWul l2005al lbT: IPapal oizou fc Ivanovl 
2005 ; llvanov fc Papaloizoull2fJfJ7l) . Inertial waves (iGreenspanl 



approach has perhaps been most successful in application 
to the circularization of binaries containing a giant star 
dVerbunt fc Phinnevlll995l ), in which case the tidal period 
exceeds the convective timescale. To study dynamical 



1968) propagate in a uniformly rotating fluid at frequencies 



ui smaller in magnitude than twice the spin frequency fi, 
as seen in a frame rotating with the fluid. The group ve- 
locity of an inertial wave is proportional to its wavelength 
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and is inclined at an angle A = arcsin |w/2fi| to the rotation 
axis. The behaviour of rays propagating within a spheri- 
cal annulus (a thick spher ical shell) is very complicated and 
sensi tive to the value of A (|Rieutord. Georgeot fc Valdettarol 
l200ll ). It generally involves the focusing of rays towards limit 
cycles known as wave attractors. Another important fea- 
ture is the existence of a critical latitude (equal to A) at 
which the rays are tangent to the core and a singularity 
is introduced into solutions of the inviscid wave equations. 
Numerical investigations indicate that the tidally forced dis- 
turbances in this frequency range are concentrated in narrow 
beams whose width diminishes as the viscosity v is reduced 
jQgilvie fc Linll2004l ). In certain intervals of uj/Q., a dissipa- 
tion rate that appears to be asymptotically independent of 
v may be achieved. However, the dissipation rate varies in a 
complicated way with the tidal frequency; presumably this 
occurs because of the sensitivity of the ray propagation to 
the value of A, and depends on the waves being reflected 
from the inner and outer boundaries. 

In an attempt to understand aspects of this behaviour, 
we have investigated a variety of reduced models. Using a 



prototyt 


fical p 


( Ogilvie 


2005) 



in a two-dimensional domain in which the rays are focused 
to wards a sing l e wav e attr actor, as occurs in the ex periments 
of iMaas et all (|l997l ) and lManders fc Maasl l|2003l ). We con- 
structed an asymptotic analytical solution and showed that 
a non-zero dissipation rate is achieved in the astrophysically 
relevant limit of small viscosity (i.e. small Ekman number). 
This limiting dissipation rate is independent of the magni- 
tude and even of the form of the small-scale damping mecha- 
nism for the waves, which is very promising for astrophysical 
applications. We confirmed the analytical solution in detail 
using numerical methods. 

The method used bv lQgilviel (|2005l ) cannot be applied 
directly to the problem of tidal forcing when the fluid occu- 
pies a spherical annulus. One technical complication is that 
the hyperbolic equation describing the structure of inviscid 
inertial waves in the meridional plane of a spherical system 
does not possess Riemann invariants, which makes it difficult 
to describe the propagation of the solution along the char- 
acteristics (which are also the ray paths described above). 
More importantly, the spherical annulus has a concave inner 
boundary with a critical latitude singularity. Although such 
a f eature could b e introduced into the problem considered 
bv lQgilvid |200i) by choosing a domain of an appropriate 
shape, separate consideration would have to be given to the 
waves generated by the singularity. 

In the present paper we consider a different idealized 
model that displays the essential features of inertial waves 
propagating within a spherical annulus and allows a de- 
tailed numerical investigation. It involves a perfectly rigid 
core surrounded by a deep ocean consisting of a homoge- 
neous incompressible fluid. Centrifugal effects are neglected, 
but the Coriolis force is considered in full, and dissipation 
occurs through viscous or frictional forces. Although the as- 
sumption of a homogeneous incompressible fluid is made for 
maximal simplicity, the model may be more or less directly 
applicable to planets or moons involving a deep ocean. We 
formulate the problem of tidal forcing in Section[2] introduce 
a further simplification in Section [3] and present numerical 



solutions in Section [4] The results are discussed in Section [5] 
and a conclusion follows in Section [6] 



2 THE TIDAL PROBLEM 
2.1 Basic equations 

We consider a simplified model of a planet in uniform rota- 
tion with angular velocity f2 = Q e z . Let (r, 6, (f>) be spheri- 
cal polar coordinates in the rotating frame. A perfectly rigid 
core occupies the region < r < aR and is surrounded by 
a homogeneous incompressible fluid with a free surface at 
r — R. We neglect the centrifugal distortion, which is of 
second order in Q. Let the mean density of the system be p 
and that of the fluid p — Pp. The surface gravity is then 



g = -nGpR. 



(1) 



We consider the long-term response of the system to 
a tidal gravitational potential Re[^(r, 8, <f>) e -1 "'], where u> 
is the (real) tidal frequency in the rotating frame. The lin- 
earized equations within the fluid arc 



2ft x u = -VW + f, 



W = ^- + $' + *, 
P 

V u = 0, 



VV = A-rGp = -A-KGp^-8{r - R), 



(2) 

(3) 

(4) 
(5) 



where u, p' , p' and $' are the Eulerian perturbations of 
velocity, density, pressure and gravitational potential, as- 
sumed to have the same form of time-dependence as the 
tidal potential, and S is the Dirac delta function. We in- 
clude a dissipative force (per unit mass) /, in the form of 
either a Navier-Stokes viscous force 



(6) 



where v is the kinematic viscosity, or a frictional damping 
force 



ft = -7", 



(7) 



where 7 is the frictional damping coefficient. The latter case 
is artificial but is mathematic ally and computationally more 
convenient. In lOgilviei (|2005r ) we found that the asymptotic 
dissipation rate in a wave attractor is identical for viscous 
and frictional forces. We set / = / v + f f in this section, 
although in the numerical analysis below we consider each 
type of dissipative force separately. 

The time-averaged energy dissipation rate is 



D = --Re 



pu* ■ fdV 



\ Re 



pWu ■ AS, 



(8) 



where S is the bounding surface of the volume V occupied 
by the fluid. The last equality follows from the linearized 
equations and the divergence theorem. For viscous and fric- 
tional damping we have 



D v = / pve i:j eij dV, 
Jv 

Dl = lJ H»| 2 dV, 



(9) 
(10) 
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where = ^(diUj + djUi) is the rate-of-strain tensor in 
Cartesian coordinates. In the viscous case, either no-slip or 
stress-free boundary conditions are required to ensure that 
there is no viscous flux of energy out of the fluid volume. 

2.2 Projection on to spherical harmonics 

To solve the linearized equations, we introduce a spheroidal- 
toroidal decompositi on and a projection on to spherical har- 
monics l|Zahnlll966ah . Thus 



a n (r)Yr(9, 



ue = r 



■E 



n[ >d8 + sm6 d<t> 
sine dd> n( ) d6 



w =^w„(r)y n m (e,0), 



(11) 

(12) 

(13) 

(14) 
(15) 



where the sums are carried out over integers n > m > 0, 
and 

1/2 



(2n + l)(n - m)\ 



T~>m t r\\ in 

P n (cos V) e 



(16) 



47r(n + m)! 

is a spherical harmonic, satisfying the orthonormality rela- 
tion 

l \Y™(6, cj>)]*Y™' (6, (ft) sin 9 d6 &</> = 5 mm i S nn i . (17) 
io Jo 

Bearing in mind the linearity of the problem, we consider 
a forcing potential in the form of a single solid spherical 
harmonic, 

V = *t(r)Yr(9,0), (18) 

for some £ > m. (Although oc r', only the value ^e(R) 
of the potential at the free surface is significant in the case 
of an incompressible fluid.) The dominant tidal potentials 
are the quadrupole components with 1 = 2, and of these the 
case m = 2 is usually the most important. The linearized 
equations couple together different values of n through the 
Coriolis force, while different values of m are decoupled. 

Projecting the equations on to spherical harmonics, we 
obtain 

(— iw + 7)o n — 2imQ.rb n 

+2flr[(n - l)<7„c n _i - (n + 2)q n +ic n+ i] 



dW n 
dr 



n(n + 1)- 



dr 



(r 2 b n ) 



(19) 



-iu> n r b n 



2imflra„ 
n(n + 1) 



+2Qr [(n - l)q n Cn-i + (n + 2)q n+1 c n+1 ] = -W n 



2a n Id/ 4 db n \ , 

h — -r- [r — - (n - l)(n + 2)b r , 

r r 2 dr V dr ' 



-iu n r c n + 2flr(q„a n -i - q n+1 a n+1 ) 

-2Qr 2 [(n - l)g„6„-i + (n + 2)q n+1 b n ^ 
"Id/ 4 dc„ 
_r 2 dr 



T r (r^)-(n-l)(n + 2) Cn 



(20) 



(21) 



H fr2 
r 2 dr 
where 



3 — (r o„) - n(n + l)b n = 0, 



0J n = LO + 



2mfi 
n(n + 1) 

2\ V2 



+ 17) 



4n 2 - 1 



qn 
n 



(22) 



(23) 



(24) 



These equa tions agree with thos e of iRieutord fc Valdettarol 
[| 19971 ) and lOrilvie fc Linl (|2004l ) in the appropriate limits. 
They are to be solved for integers n > m, but in practice 
the system must be truncated. 

The potential perturbation satisfies 



l__d_ 

r 2 dr 



dr 



n( - n + V tf n = -^Gp^5(r - R), (25) 
r 2 vjj 



and the solution that is well behaved as r — » and as r — *• 00 
is 

, A-KGpRa n {R) f ( r /R) n , < r < R, 

$ " - 2n + l X i (r/R)-^\ r > R. (26) 



2.3 Boundary conditions 

Since the inner boundary is impermeable, a„ = at r = aR. 
At the outer boundary, the normal stress at the perturbed 
free surface vanishes. Thus —p' n —pg(a„/iu})+2pv(da n /dr) = 
0, i.e. 



Wn - 2w 



da n 
dr 



+ 



3/3 
2n + 



1UJ 



(27) 



at r = R. In the limits of low frequency and low viscosity, 
lu 2 <C gl/R and v <C gR/u, the first two terms are negligible, 
and this equation amounts to specifying the radial velocity 
at the surface. Indeed, it amounts to saying that the radial 
displacement of the surface is equal to the simple equilibrium 
tide —ty/g, with a correction (the term involving 0) due to 
the self-gravity of the fluid. 

When viscosity is included, additional boundary condi- 
tions are required. We assume that the tangential stresses 
vanish at the inner and outer boundaries, so 



dc n 



an db n _ 

r 2 dr dr 







(28) 



at r = aR and r = R. The inner boundary could alterna- 
tively be treated using the no-slip conditions b n — c n — 0, 
which would result in a prominent Ekman or oscillatory 
boundary layer. 



2.4 Dissipation 

The time-averaged energy equation for this system has the 
form 



dF A W 

dr 

where 



(29) 



F=lr 2 pRe^a* n (w n - 2u^ 

n 

—n(n + l)u ^b n (a„ 



dc„ 



+ r — — J + r c n 

dr J dr 



(30) 
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is the radial energy flux and 
d v — \-r 2 ' pv S ' n(n + 1 



r dr 



. dCn 2 

dr 



I dr I 



d i = '^ r2 / 3 7^ |a«| 2 + n{n + l)r 2 (\b n \ 2 + \c n \ 



(32) 



are the viscous and frictional dissipation rates per unit ra- 
dius. The integrated energy equation, together with the 
boundary conditions, relates the total dissipation rate to 
the surface forcing: 



D 



D v +D f 

l-R 

(d v + df) dr 

aR 

l -pR 2 Re[a e {R)*q e {R)}. 



2.5 Frictional non-rotating problem 



(33) 
(34) 



In the absence of rotation the different spherical harmonics 
are decoupled. When only a frictional force is present, the 
solution satisfying the inner boundary condition is 



at = A 



A 
r 



21+ 



e+2' 



+ 



'(*) 

y 2e+1 /R\ 

e+i \r) 



e+2' 



We — (iw — ^)Ar 



1 / r y- 1 a 



+ 1 \r) 



(35) 
(36) 
(37) 



where A is a complex coefficient, and all other components 
vanish. The outer boundary condition determines the value 
of A according to 

fl a 2i+1 ' 



+ 



i - 



3/3 



,2£ + 1 
and the dissipation rate is 



2- A (1-a 
ltd y 



2l+\\ 



(38) 



D 



2£+l\ 



£+1 



(39) 



We specialize to the case of a quadrupolar tide (I = 
2). The dissipation rate can then be converted into a tidal 
quality factor Q using the relation 

15 R\uj\\^2{R)\ 2 



D 



8Q' 



2-kG 



(40) 



(Here Q' = 3Q/2fc2 is the modified tidal quality factor, k2 
being the second-order potential Love number. For a homo- 
geneous fluid body &2 = 3/2 and Q' — Q.) The solution of 
equation (|38|) is 



A 



J^2(R) 

9 



B, 



(41) 



where B is given by 



B-i = (1 _ a «)(i _ "fl _ (1 + |a 5 )^^l (42) 

and embodies various corrections (due to the self-gravity of 
the fluid and the non-zero tidal frequency) to the simple 
equilibrium tide. Then 

\B[ 2 
5 



(43) 



Note that a resonance occurs with the I = 2 surface gravity 
mode (f mode) as the tidal frequency is increased into the 
vicinity of the dynamical frequency (g/R) 1 ^ 2 ; in the coreless 
case (a = 0, (3 = 1) this occurs at uj 2 = Ag/5R. 



2.6 Viscous non-rotating problem 

When the fluid is non-rotating and has only viscous dissipa- 
tion, the general solution of the linearized equations in the 
case £ — 2 is 

a2 = 4 [ AlT ' 5 + A 2 + M 3 - 3ifcr - fcV) e ifer 



+Ai(3 + 3ifer - fcV) e -1 * 



(44) 



-g [3Air 5 - 2A 2 - A 3 (6 - 6ifcr - 3fcV + ifcV) e ifc 



-A 4 (6 + 6ifcr - 3fcV - ifcV) e" ifcrl 



W 2 = ^(3Air°-2A 2 ), 
where Ai, A2, A3 and A4 are complex coefficients and 



(45) 
(46) 



k ■■ 



(?) 



1/2 



(47) 



Note that the terms involving an exponential correspond to 
viscous boundary-layer solutions, wh ile the other te rms are 
solutions of the inviscid problem (cf. lRieutordlll99~lh . 

With stress-free inner and outer boundary conditions, 
the solution for the velocity in the limit of small u is ev- 
erywhere close to the inviscid one. The boundary-layer solu- 
tions, which decay rapidly with distance, intervene near the 
boundaries to adjust the shear rate to match the stress-free 
conditions. Viscous dissipation is dominated by the region 
away from the boundaries, for which the shear rate is essen- 
tially that of the inviscid solution. We then find 



D « -(3 + 5a 3 - 8a 10 )pisR\A\ 2 , 
6 



(48) 



with A as given by equations (141|l and (|42|l (with 7 = 0). 
Thus 



Q' 



2\B\ 2 (1 + fa 3 



§a 10 )/3^ 



gR 



(49) 



With a no-slip inner boundary condition, viscous dissi- 
pation would be dominated by the shear layer on the inner 
boundary. In this case the dissipation rate would be approx- 
imately proportional to v 1 ^ 2 , rather than v, in the limit of 
small viscosity. 



2.7 Tidally forced inertial waves 

When rotation is included, the linearized equations support 
inertial waves for frequencies in the range —2 < u/Sl < 2. 
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This can be seen by setting the dissipative force / to zero 
and deducing the Poincare equation 



(50) 



to 2 y 2 W - (2SI • V) 2 W = 

for the modified pressure perturbation l|Greenspanl fl968). 
When rendered in cylindrical or spherical polar coordinates 
for solutions with the azimuthal dependence e 1 " 1 '*, this be- 
comes a second-order partial differential equation that is of 
hyperbolic type for < lu 2 < 4Q 2 . The characteristics are 
straight lines in the meridional plane inclined at an angle 
A = arcsin |cj/2fl| to the rotation axis. 

Although the Poincare equation is homogeneous, the 
waves are forced through the inhomogeneous outer bound- 
ary condition. Since a hyperbolic equation is generally ill 
posed when supplied with boundary conditions on a closed 
surface, the inviscid problem generally has no solution. Ei- 
ther viscosity or friction formally removes the hyperbolic 
character of the equation and leads to a well posed problem. 
The initial-value problem for an inviscid fluid is also well 
posed. 



3 RADIALLY FORCED OSCILLATIONS 
3.1 Radial forcing 

We have seen that, in the limits of low frequency and low 
viscosity, the outer boundary condition amounts to specify- 
ing the radial velocity at the surface. We therefore consider 
a slightly simpler problem in which the free outer surface 
is replaced by a boundary (having zero tangential viscous 
stress) on which the radial velocity u r is prescribed and is 
in the form of a single spherical harmonic, i.e. 



U5 n t 



(51) 



at T = R, where U is an arbitrary (real) amplitude. The 
dissipation rate is then 

D = -^pR 2 URe[W e (R)] . 



3.2 Frictional problem 



(52) 



When only a frictional force is present, we can eliminate 
variables in favour of the radial velocity to obtain 



( n- 1 \ gj + fn + 2\ 
\ n J w„-i \n + 1 J 



~2 

1n + l 



+ 



+ 



Wn-i Vn + 1/ LOn+i 4n(n + 1)Q, 2 
x(d y - n + l)(d y + n + 2)a n 

q n q n -\ ^ — n + l)(d y — n + 3)a„_ 2 
Wn-l 

UJn+1 



(d y +n + 2)(d y +n + 4)a n+2 = 0, 



(53) 



where y = ln(r/_R). 

We truncate the system at spherical harmonic degree L 
by setting a„ = for n > L. The truncated problem consists 
of a finite system of homogeneous linear ordinary differential 
equations (ODEs) with constant coefficients. This possesses 
a basis of exponential solutions with a n oc exp[(fc — l)y] oc 
r fe_1 for suitable values of k, in which the values of a n are 
related by a three-term recurrence relation of the form 

X n a n -2 + Y n a n + Z n a n+2 = 0. (54) 



In the case £ = m = 2, for example, the permitted velocity 
components are n = 2, 4, 6, . . . , L and the coefficients of the 
recurrence relation satisfy X2 = Zl = 0. The system of 
ODEs is of order L and its exponential solutions are of the 
following form: 

(i) Solutions with k = 2, 4, 6, . . . , L. These can be nor- 
malized by setting a 2 = 1, and iteration of the recurrence 
relation yields 04, ae, ■ ■ ■ , o-h- Beyond this point, = 
because Xk = = 0. (These are exact solutions of the 
untruncated problem.) 

(ii) Solutions with k = — L — 1, —L + 1, — L + 3, . . . , —3. 
With the normalization a_L = 1, the recurrence relation 
yields ol-2, ol-4, • ■ • , d-fc-i- Then a_k-3 = because 
V_fc_i = = 0. (These solutions are exact only in 
the truncated problem.) 

In a related approach, iRieutordl fl991) obtained analytical 
solutions of the viscous problem in terms of Bessel functions 
and polynomials in r. 

It is possible to solve the problem of forced inertial 
waves by writing the solution as a linear combination of 
the exponential solutions and applying the boundary condi- 
tions to determine the unknown coefficients. However, the 
solutions with large values of \k\ vary by many orders of 
magnitude between the inner and outer boundaries, except 
possibly in the case of a thin shell, and this gives rise to 
a highly ill conditioned problem. It would be interesting to 
investigate whether this problem can nevertheless be solved 
numerically for large values of L. In this paper, however, 
we mainly solve the ODEs using the Chebyshev collocation 
method, as discussed in Section f4. II below. 



3.3 Normal modes of the inviscid problem 

Do there exist normal modes for inertial waves in a spheri- 
cal annulus? Since we are now considering a radially forced 
problem, normal modes would be solutions of the same prob- 
lem but with the homogeneous outer boundary condition 
a n — 0. These are free oscillations within rigid boundaries, 
and should agree with the oscillation modes of a fluid with a 
free outer surface in the limit of frequencies much lower than 
those of surface gravity waves. We consider a fluid without 
dissipative forces; if it possesses a normal mode then the 
frictional problem admits the same mode but with an addi- 
tional damping rate of 7. (Viscous fluids certainly do possess 
normal m odes, but they tend to be loc alized around wave 
attractors; iRieutord &: Valdettard[l997l .') 

There do exist special, purely toroidal modes, for which 

W = -^T' ^ 

Cm = r m_1 , (56) 

W m+1 = ~2m{lq m+1 r m+1 , (57) 

while all other components vanish l)Rieutord fc Valdettard 
Il997t) . These solutions exist independent of the radial extent 
of the annulus. Note that there is only one such mode for 
each azimuthal wavenumber m. 

Normal modes involving some radial motion can be 
sought by expressing them in the basis of exponential solu- 
tions described above. Setting the determinant of the matrix 
of coefficients to zero yields (in the case m = 2) a polynomial 



6 Gordon I. Ogilvie 



equation of degree L(L — 2)/2 in oj/fi. Numerical investiga- 
tion suggests that the roots are all real and lie in the interval 
(—2,2). These solutions presumably give rise to a complete 
set of normal modes for the truncated problem. However, 
the frequencies and eigenfunctions of the modes do not ap- 
pear to converge as L is increased, which suggests that the 
spherical annulus does not possess normal modes other than 
the purely toroidal ones. 

In the speci al case of a corel ess planet (a = 0) normal 
modes do exist l|Greenspanl 1 19681 ) . The simplest symmetric 
modes for m — 2, for example, are 



»(r 2 -R 2 )5 n2 , 



fi 




-1.2319, 0.2319. 



(58) 



(59) 



It might be though that such modes could be excited reso- 
nantly by tidal forcing. However, the exact solution of the 
radially forced problem for £ = m — 2 in a coreless planet is 



a n = — 8n2 



(60) 



and does not allow for any such resonances. (In cylindri- 
cal polar coordinates (s, 4>,z) this solution is u s — s/R, 
— is/ R, u z = 0.) This is also a solution of the nonlinear 
problem, because u ■ Vu = 0. In a similar way, the Roche- 
Riemann ellipsoids provide the nonlinear solutions of certain 
problems of tidal forcing in the case of a homogeneous in- 
comp ressible planet without a core (|Goodman fc Lacknerl 
12003 ). 



4 NUMERICAL ANALYSIS 
4.1 Method of solution 

We focus on the i — m — 2 type of quadrupolar forcing that 
is of greatest importance in applications. For the viscous 
problem, we solve equations (|19p ~ (|22fl and the appropriate 
boundary conditions by t he ps eudospectral method as de- 
scribed in lOgilvie fc Linl (|2004h . This makes use of Cheby- 
shev collocation on a grid of N + 1 points in r (the Gauss- 
Lobatto nodes) . The system of equations is also truncated at 
spherical harmonic degree L. The problem involves a large 
block-tridiagonal matrix and the solution is obtained by a 
standard direct method. The total dissipation rate is com- 
puted by two independent methods: from the direct dissi- 
pation integral (equation [33} using a Chebyshev quadrature 
formula, and from the rate of working at the surface (equa- 
tion [34] or [52] as appropriate). 

For the frictional problem, we apply Chebyshev colloca- 
tion to equation (|53p and the appropriate boundary condi- 
tions. This has the advantage that the blocks of the tridiago- 
nal matrix are smaller because the other variables have been 
eliminated in favour of the radial velocity, and therefore a 
higher resolution can be reached. (In this case the Cheby- 
shev variable is taken to be linearly related to y rather than 
r.) 

To validate these methods, we checked that the solu- 
tions found by the exponential technique described in Sec- 
tion !3.2l which works very well for modest values of L, agree 
with those obtained by the pseudospectral method. We ver- 
ified numerically the analytical expressions obtained in Sec- 



tions [23] and [23] for the dissipation rate in the non-rotating 
tidal problem. We also confirmed that the dissipation rates 
calculated by the integral and surface formulae agree well in 
all cases. 



4.2 Overview of results 

In Figs \T\ and [5] we present an overview of the behaviour 
of the dissipation rate with tidal frequency in the radially 
forced problem with a frictional force, for different values 
of the fractional core radius a and the frictional damping 
coefficient 7. These results were obtained with numerical 
resolutions of L = N = 200 or 300, which is adequate for this 
kind of overview. A higher resolution is generally required 
for smaller values of 7. 

Outside the range —2 < w/fi < 2 the dissipation rate 
is essentially independent of u> and proportional to 7. No 
inertial waves are excited and the response is smooth and 
uninteresting. The dissipation rate is greater in a thin shell 
because of the larger value of the parameter B. 

Within the range — 2 < tj/fi < 2 the dissipation rate is 
en hanced through the e xcitation of inertial waves, as found 
bv lQgilvie fc Linl <|2004h . There is a complicated dependence 
on frequency. In the frequency intervals in which the dissi- 
pation rate is most enhanced, it appears to be roughly in- 
dependent of the frictional damping coefficient. 

This enhancement of the dissipation rate increases sys- 
tematically with the size of the core. In the case of a very 
thin shell (a — 0.99), we observe instead a strong resonance 
with a Rossby wave or r mode at ui = —fi/3, as described 
below. 



4.3 Resonance with a Rossby wave or r mode in a 
thin shell 

The limit of a thin shell, 1 — a = e -C 1, is of some in- 
tere st and can b e related to problems studied by Laplace 
(see lLamblll932l . and references therein). In the absence of 
viscosity, motions are possible in which, to a first approx- 
imation, a n varies linearly with r, vanishing at the inner 
radius, and is much smaller than rb n or rc n . Under these 
circumstances the 'traditional approximation' and Laplace's 
tidal equations are applicable. In this limit the shell pos- 
sesses free modes of oscillation which depend on the dimen- 
sionless para meter 4fi 2 R 2 /gh, wher e h = eR is the depth 
of the ocean |Longuet-Hi ggins 1968:). When this parameter 
is much less than unity, simple analytical solutions emerge 
with frequency lj = — 2mfi/n(n + 1) for integers n > m. 
These are the planetary or Rossby waves, which are well sep- 
arated from surface gravity waves in the limit under consid- 
eration. Related solutions can be identified in stably strati- 
fied stellar atmospheres, wh ere they are known as r modes 
jPapaloizou fc Pringlelll978j ). 

The resonance described above occurs with the (m = 
2, n = 3) Rossby wave or r mode at U) = —fi/3. (Note that 
this is not equivalent to the special, purely toroidal mode dis- 
cussed in Section [3. 31 which has a frequency of lo = — 2fi/3 
and involves no radial motion.) The appropriate balance in 
equations (|21|l and (1221) is then 

- \LJ n c n - 2fi[(n - l)q n b n -i + (n + 2)q„+ib n+1 ] W 0, (61) 



Tidal dissipation in rotating fluid bodies 7 









o 





-3-2-10 1 2 3 

Figure 1. Variation of the total dissipation rate (in units of pR 3 U 2 fi) with tidal frequency, for a fluid with a frictional force. The radial 
velocity is forced at the outer boundary in the form of an I = m = 2 spherical harmonic of amplitude U. The fractional core radius is 
a = 0.1, 0.2, 0.3, 0.4 and 0.5 (panels from top to bottom). The frictional damping coefficient 7 is given by 7/n = 10~ 2 , 10~ 2 - 5 and 10~ 3 
(red, green and blue curves, from top to bottom for the most part). In all figures, 'log' denotes a base-10 logarithm. 
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Figure 2. Continuation of Fig.[T]for fractional core radii a = 0.6, 0.7, 0.8, 0.9 and 0.99 (panels from top to bottom). Note the different 
vertical scale in the last case. 
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Figure 3. Comparison of tidal and radially forced problems. The bottom (blue) curve shows the total dissipation rate (in units of 
pR 3 U 2 n) for a = 0.5 and j/Q, = lO" 3 , as in the bottom panel of Fig. [T] The middle (green) and top (red) curves show the dissipation 
rate for the tidally forced problem in a planet with dynamical frequency (g/R) 1 / 2 = 5fl and 3 £7, respectively, and with mean density 
equal to the fluid density. The latter curves are vertically offset by 2 and 4 units, respectively, for clarity. The comparison is made by 
equating the forced radial velocity amplitude U with u)(5/2)('$> /g). 



n(n +!)&„ ~ 0. 



(62) 



dr 

The boundary conditions determine that da n /dr = 
(U/eR)S n £ and we easily obtain b n and c n from the above 
equations. In the case I = m = 2 investigated here, the 
toroidal velocity component C3 is excited and becomes large 
when U13 is small, which occurs close to u)/Q, = —1/3. The 
resulting dissipation rate is 



D 



pjR 3 U 2 
12e 



1 + 



32fi 2 
63|ci>3 



(63) 



This expression provides an excellent fit to the numeri- 
cally determined dissipation rate in the frictional problem 
for a — 0.99. Note that D is inversely proportional to the 
thickness of the shell. This result occurs because the fast 
barotropic flow (i.e. horizontal motion independent of depth) 
excited in a thin shell is efficiently damped by the artificial 
frictional force in this model; a different behaviour would oc- 
cur in a viscous fluid. Nevertheless, it is well known from the 
case of the Earth that highly efficient tidal dissipation can 
occur in a shallow ocean. The case of the Earth's ocean is 
of course enormously complicated by its irregular shape and 
depth. Mode conversion and turbulence provide channels for 
dissipation. 

4.4 Comparison of tidal and radially forced 
problems 

The original tidal problem, described in Section gives re- 
sults very similar to those of the simplified, radially forced 
problem. Some additional parameters are required to set 
up the tidal problem: the ratio of the dynamical frequency 
(g/R) 1/2 to the spin frequency f2, and the ratio f3 of the 
fluid density to the mean density of the planet. In Fig. [3] 



we compare the dissipation rates in the tidal and radially 
forced problems, with a — 0.5, (3 = 1 and = 10 . In 
order to make the comparison we assume that the radial ve- 
locity at the surface is determined by equation (|27[) in the 
low-frequency limit in which the W„ term (and also the vis- 
cous term) is neglected. This is equivalent to saying that the 
radial displacement of the surface is — (5/2)^ /g, where the 
factor of 5/2 comes from the self- gravity of the fluid. 

When the planet is slowly rotating, in the sense that 
Q -C (g/i?) 1//2 , the agreement is excellent and shows that 
the excitation of inertial waves is not affected by the freedom 
of the outer boundary. For more rapidly rotating planets the 
influence of surface gravity waves becomes apparent in this 
range of frequencies. For such rapidly rotating planets the 
centrifugal distortion of the fluid ought to be taken into 
account. 



4.5 Investigation of a restricted frequency interval 

In order to investigate the response in greater detail, we 
restrict our attention to the case of a fractional core radius 
a = 0.5 and to a certain range of frequencies, 1 < uj/Q < 1.2. 
This interval is chosen partly because, as described below, it 
contains two relatively simple wave attractors. Fig. [4] shows 
an expanded view of the bottom panel of Fig. [T] using a 
higher frequency resolution. A higher numerical resolution 
(up to L = N = 800) is also used in order to ensure that 
the results are adequately converged for the appearance of 
this figure. 

Fig. [5] shows equivalent results obtained for a viscous 
fluid. In this case a smaller numerical resolution (up to 
L = N = 400) is achieved and the results are adequately 
converged except possibly in some cases for the lowest vis- 
cosities. 
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Figure 4. Expanded view of the bottom panel of Fig. [T] i.e. for 
a = 0.5 and y/Cl = 10~ 2 , 10~ 2 ' 5 and 10" 3 . 








-5 



1.00 1.05 1.10 1.15 1.20 

Figure 5. Equivalent of Fig. [4]for a viscous fluid with v/R 2 Q = 
10 -4 , 10 — 5 , 10 -6 , 10 — 7 and 10 -8 (red, green, blue, magenta and 
cyan curves). 

We also show the variation of the dissipation rate with 
the coefficient of friction or viscosity for three frequencies in 
this interval, lj/Q. = 1.05, 1.10 and 1.15, in Figs E] and 
To obtain adequately converged results for the lowest values 
of 7 or v, a very high resolution (up to L = N = 1600) is 
required. 

Perhaps the simplest behaviour occurs close to cj/fi — 
1.10. Here the dissipation rate converges to a non-negligible 
value as either 7 — -> or v — > 0. The same limiting value is 
obtained in the frictional and viscous problems and depends 
smoothly on to within a certain interval. This behaviour is 








— 4 r .... 1 .... 1 .... 1 .... 1 .... 1 .... • 

-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 
log( 7 /Q) 

Figure 6. Frictional dissipation versus frictional damping coeffi- 
cient for a = 0.5 and lj/Q = 1.05, 1.10 and 1.15 (red, green and 
blue curves, from top to bottom). 
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Figure 7. Viscous dissipation versus viscosity for a = 0.5 and 
lj/H = 1.05, 1.10 and 1.15 (red, green and blue curves, from top 
to bottom). 



closely reminisc ent of that fou nd in the analysis of a simple 
wave attractor We will see below that two 

distinct types of wave attractor are active in this range of 
frequencies. 

The maximum dissipation rate occurs close to uj/Q. — 
1.05. In the frictional problem D does not appear to converge 
as 7 — > 0, although the ultimate behaviour is not clear from 
the limited range of parameters that can be investigated 
numerically. In the viscous problem the dissipation rate gives 
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Figure 8. Structure of the velocity perturbation \u\ in a merid- 
ional quarter-plane for the case a = 0.5, uj/Q = 1.05, 7/H = 
10 -3 , L = N = 1600. A linear colour table is used. 

the appearance of convergence as v — > but it may decrease 
below v/R 2 n = 10" 9 . 

A much lower dissipation rate is found close to ui/Q = 
1.15. Here D appears to converge to a very small value as 
v — > 0. It may do so in the frictional problem at lower values 
of 7 than are achieved here. 

4.6 Spatial structure of the response 

We now examine the spatial structure of the response in the 
meridional plane. FigslB nTOl show the magnitude of the veloc- 
ity perturbation, \u\, for frequencies U)/Cl = 1.05, 1.10 and 
1.15, in the frictional problem with 7/fi = 10~ 3 . In each case 
the disturbance is concentrated near selected inertial-wave 
rays, which are characteristics of the Poincare equation. 

Relatively simple wave attractors exist in this range 
of frequencies. The symmetrical pair of attractors involv- 
ing four reflections on the outer boundary, illustrated in the 
upper part of Fig. [TTI exists for 1.076 < u)/Q < 1.148; a sec- 
ond pair of attractors, involving six outer reflections, exists 
for 1.091 < ui/Q < 1.165. The focusing power of an at- 
tractor depends on frequency in a systematic way within 
the bandwidth of the attractor, as its shape is adapted 
jRieutord, Georgeot fc Valdettarol l200lf) . The first one is 
strongest at the right-hand end of its bandwidth, while the 
second is strongest at the left-hand end. 

The response at frequency 1.10 (Fig. [9| shows that the 
two types of attractor that exist at this frequency are both 
activated. The wave energy is concentrated around the at- 
tractors in a region whose width is proportional to 7. This 
behaviour, and the convergence of the dissipation rate as 
7 — » or v — > 0, are entirely analogous to that found by 
lOgilviei (|2005l ). 

Also seen in this case, however, is a ray emerging from 
the critical latitude on the inner boundary. This ray is indeed 
the dominant feature apparent at frequency 1.15 (Fig. I10|) 




Figure 9. As Fig.[8]but for uj/Q = 1.10. 




Figure 10. As Fig.|8]but for u/Q = 1.15. 



where the dissipation is relatively weak. The attractor with 
six outer reflections is present in the solution at frequency 
1.15 but is weakly focusing and is not powerfully activated. 
The dissipation rate in the viscous problem appears to con- 
verge to a very low value as v —* at this frequency. 

The behaviour at frequency 1.05 (Fig. |SJ), where the 
dissipation is strongest, is harder to explain. The ray dy- 
namics in the vicinity of this frequency is very complicated. 
Although extremely long attractors may exist (e.g. one at 
frequency 1.05 involving 584 outer reflections) they are of 
no practical significance. However, there is a tendency for 
rays to be concentrated temporarily into the region between 
the equator of the inner boundary and the critical latitude 
on the outer boundary, where the wave energy is seen to 
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Figure 11. Wave attractors in a spherical annulus with a = 0.5 
at frequency uj/CI = 1.10. 



be prominent. The rays escape from this region and reenter 
it repeatedly. It is difficult to identify features of the ray 
dynamics that lead to a preference for the frequency 1.05. 

The prominence of rays emerging from the critical lat- 
itude on the inner boundary has been noticed before (e.g. 



ltuae o n tne inner Boun dary nas 
lTilgnerlll999l ; lOgilvie fc Lhj|2004h 



4.7 Frequency-averaged dissipation rate 

There is a strong systematic dependence of the dissipation 
rate on the size of the core. This effect is examined in Fig. 1121 
where we plot an average dissipation rate versus a. The 
average is obtained from the data for Figs[T]and[2] by taking 
the arithmetic mean of the dissipation rates evaluated at 
800 equally spaced frequencies in the range —2 < uj/Q. < 2. 
The 'baseline' dissipation rate, which occurs for frequencies 
outside this range, is subtracted from the average, in order 
to isolate the contribution due to waves. There may be a 
systematic error in this procedure owing to the inadequate 
sampling of narrow peaks in the curve. 

The resulting values are approximately independent of 
7 and increase approximately as the fifth power of the core 
radius. Interesti ngly, this dependence is pre dicted by the 
recent analysis of lGoodman fc Lacknerl l|2008h , who consider 
the scattering of the equilibrium tide off the core, although 
they do not predict the intricate variation with frequency. 
A departure from this behaviour occurs as a approaches 1, 
owing to the toroidal-mode resonance in a thin shell. 



5 DISCUSSION 

5.1 Interpretation of enhanced dissipation 

The dynamics of inertial waves in a spherical annulus is 
very rich a nd a great de al remains to be explained. In ear- 
lier work (|Ogilviel [iooHh we considered a simpler type of 



Figure 12. Estimates of the frequency-averaged dissipation rate 
over the interval —2 < w/Q < 2, minus the 'baseline' dissipation 
outside this range, versus fractional core radius, for the frictional 
problem with 10~ 3 < y/Cl < 10~ 2 . The dotted line represents 
0.11a 5 . 



domain in which all wave energy is focused towards a sin- 
gle wave attractor and a dissipation rate is obtained that is 
asymptotically independent of the viscosity or other small- 
scale damping mechanism for the waves. In a spherical an- 
nulus, the closest parallel to this behaviour occurs in sit- 
uations such as that illustrated in Fig. [§] although in fact 
more than one attractor is activated there. The spherical 
annulus is more complicated because simple attractors oc- 
cupy only a limited range of frequencies, and also because 
of the existence of a critical latitude at which rays are tan- 
gent to the inner boundary. In fact these features are related 
to each other, because the range of frequencies in which an 
attractor exists is bounded by situations in which its rays 
appro ach the critical latitude on the inner and outer bound- 
aries l|Rieutord. Georgeot fc Valdettaroll200ll 'h 

Both attractors and critical latitudes concentrate wave 
energy and promote enhanced dissipation. A smooth beam 
of waves reflecting from the inner boundary generates a 
singularity in an inviscid fluid, as illustrated in Fig. 1131 
This effect depends smoothly on frequency in the range 
—2 < uj/Q < 2. When the singular beam emanating from the 
critical latitude is resolved by viscosity, enhanced dissipation 
occurs but this alone does not account for the behaviour seen 
close to uj/Q. — 1.05 in Fig. [5] it is more suitable as an ex- 
planation of the behaviour close to u/fl = 1.00, where D 
scales approximately with u 1 ' 2 . 

The behaviour of rays alone is not sufficient to explain 
the shape of the dissipation curve. This can be seen from the 
fact that D depends on the sign of ui/Sl and on the value of 
m (not shown here), whereas the rays are insensitive to these 
properties. Even in the case of a simple wave attractor, the 
full calculation of the asymptotic dissipation rate depends 
on a detailed a nalysis of the forcing accumulated around 
each ray circuit l|Ogilvieil2005h . 
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Figure 13. Formation of a singularity at the critical latitude 
through the reflection of rays (equally spaced lines incident from 
the top right) from the inner boundary. 

Some time after the calculations described in this 
paper were completed, w e received a preprint from 
iGoodman fc Lacknerl (|2008l ) . They treat the equilibrium 
tide as a large-scale inertial wave that reflects from the in- 
ner boundary. The intense beam emanating from the neigh- 
bourhood of the critical latitude produces enhanced dis- 
sipation, occurring in their analysis through nonlinearity. 
This theory has an attractive simplicity to it and appears 
to predict the correct scaling of the frequency-averaged dis- 
sipation rate on the size of the core. It does not predict 
a complicated frequency-dependence because further reflec- 
tions are not considered. Indeed, as is mentioned below, the 
various physical circumstances in planets and stars may not 
always permit the multiple reflections that occur in the sim- 
ple model considered in the present paper. 

After this paper was submitted, we also received a 
manuscript from M. Rieutord and L. Valdettaro, who have 
considered a similar problem and have been able to relate 
some aspects of the frequency-dependence of the dissipation 
rate to the properties of wave attractors, viscous normal 
modes and critical latitudes. 

5.2 Astrophysical and planetary applications 

The dissipation rates D (expressed in units of pR 3 U 2 Q.) plot- 
ted in various figures in this paper can be converted into 
modified tidal quality factors Q' according to the relation 

Q' 5 Pl 5 ^ g/R \pR^U 2 9. ) 

For a fixed value of w/f2 in the range of inertial waves, the 
efficiency of tidal dissipation through this mechanism there- 
fore scales with the square of the ratio of the spin frequency 
to th e dynamical frequency of the body (|Ogilvie fc Lin! 
|2004| ). This fact should always be borne in mind when com- 
paring solar system planets with hot Jupiters and hot Nep- 
tunes, which are expected to rotate more slowly as a result of 
tidal evolution. The wide variation in D indicates that a very 



broad range of Q' values can occur through this mechanism. 
The strong frequency-dependence of Q' can have important 
consequences for tidal evolution, for example in shaping the 
satellite systems of the giant planets. 

We note that the tidal quality factor provides 
only a conventional and, for some purposes, convenient 
parametrization of the efficiency of tidal dissipation, which 
can involve complicated processes of wave excitation, prop- 
agation and damping. It depends as much on the capacity of 
the system to respond to oscillatory forcing as on its ability 
to dissipate the response. 

5.3 Topics for further investigation 

The problem studied in this paper is intentionally simplified 
in order to try to isolate an already complicated physical 
process. We mention here some of the topics that require 
further investigation. 

The quality of wave reflection from the boundaries of 
the fluid region supporting inertial waves should be anal- 
ysed. Depending on the type of planet or star being consid- 
ered, these boundaries may be free surfaces, interfaces with 
solid or denser fluid material, or borders between convective 
and radiative zones. Internal reflection of inertial waves from 
discontinuities such as a putative first-order phase transition 
in giant planets are also potentially of importance. 

The propagation of inertial waves could be affected by 
differential rotation, magnetic fields, buoyancy effects and 
the centrifugal distortion of the body. It would also be of 
interest to study an initial-value problem rather than as- 
suming a response that is harmonic in time. 

Dissipation mechanisms for inertial waves at small 
scales include viscous damping, interaction with turbulent 
convection, Ohmic dissipation through coupling with mag- 
netic fields, and nonlinear dissipation through wave break- 
ing or parametric instability. While the analysis of a sim- 
ple wave attractor indicated a total dissipation rate that is 
asymptotically independent of the small-scale processes, a 
more complicated picture arises from the present paper and 
it may be necessary to understand the damping mechanisms 
in greater detail. 

The dependence of the dissipation rate on the size of the 
core should be investigated in more realistic models where 
the fluid is compressible and the core is deformable. The 
role of resonances with global normal modes (if they exist) 
in compressible models needs to be clarified. The relevance 
of boundary layers on the surface of the core, obviated here 
through the use of stress-free boundary conditions, is wor- 
thy of consideration. Dissipation within the core itself is yet 
another area of uncertainty. 

5.4 Elliptical instability 

A quite different paradigm for tidal dissipation in rotating 
planets and stars co uld emerge thro ugh a consideration of el- 
liptical instabilities (Kcrswcll 2002, and references therein). 
If the linear tidal response is considered to provide merely 
an elliptical distortion of the streamlines of the rotating 
fluid, as is true in the case of a full sphere of incompress- 
ible fluid, then the secondary instabilities of the elliptical 
flow can produce turbulence and enhanced tidal dissipation. 
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Analytical, numerical and experimental investigations have 
been moti vated mainly by the case of the Earth's core but 
also by Io (Kerswell & Malkus 199 a) and by accretion discs 
in binary stars ([Good man 19931; iLubow. Pringle fe Kerswelj 
1 19931 ; iRvu fe Goodmanl I 19941 ). This nonlinear mechanism 
may be more relevant for tides in extrasolar planets, while 
the much weaker tides in solar-system planets may depend 
on linear mechanisms such as those described in this paper. 
Since the linear response in a spherical annulus departs sig- 
nificantly from a simple elliptical distortion, it would be of 
interest to understand the nonlinear outcome in the presence 
of a core. 



6 CONCLUSION 

In this paper we have studied the tidal forcing, propaga- 
tion and dissipation of linear inertial waves in a rotating 
fluid body. The intentionally simplified model involves a per- 
fectly rigid core surrounded by a deep ocean consisting of 
a homogeneous incompressible fluid. Centrifugal effects are 
neglected, but the Coriolis force is considered in full, and dis- 
sipation occurs through viscous or frictional forces. We in- 
troduced a further simplification by replacing the free outer 
surface with a boundary on which the radial velocity is spec- 
ified, which closely mimics the effects of tidal forcing at low 
frequencies. Various analytical results provide tests of the 
numerical methods for calculating forced linear waves. 

The dissipation rate exhibits a complicated dependence 
on the tidal frequency and generally increases strongly with 
the size of the core. In certain intervals of frequency, effi- 
cient dissipation is found to occur even for very small values 
of the coefficient of viscosity or friction, which is promising 
for astrophysical and planetary applications. In restricted 
intervals, the wave energy is focused towards relatively sim- 
ple wave attractors and a well defined asymptotic dissipa- 
tion rate is achieved. However, a more typical behaviour is 
that the inertial waves propagate in a very complicated way 
around the fluid annulus. The critical latitude on the inner 
boundary plays an important role in the solutions. While 
the dissipation rate is enhanced in a strongly frequency- 
dependent manner, it may not converge in the limit of small 
viscosity in the same way as for a wave attractor. In the 
limit of a thin fluid shell, the dissipation rate can be greatly 
enhanced through a traditional type of resonance with a 
toroidal mode (r mode or Rossby wave). 

The model adopted in this paper is deliberately oversim- 
plified, although it may be more or less applicable to planets 
or moons involving a deep ocean. We have pointed out nu- 
merous avenues for further investigation, including aspects 
of internal wave propagation, reflection and dissipation, and 
nonlinear behaviour such as the elliptical instability. 
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